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ABSTRACT 


Constraint  has  been  an  important  consideration  in  fracture  mechanics  from  the  earliest 
work  that  was  done  to  develop  the  1974  version  of  the  ASTM  Standard  E399.  Stringent 
thickness  and  size  requirements  were  placed  on  the  test  samples,  in  terms  of  the  measured  plastic 
zone  size.  These  requirements  often  invalidate  the  results  in  practical  application  sizes.  Similar 
size  requirements  have  been  incorporated  into  more  recent  elastic-plastic  fracture  test  standards 
like  ASTM  E813  and  ASTM  El  1S2,  even  though  experimental  data  has  not  shown  a  clear  and 
direct  relationship  between  specimen  size  and  the  resulting  initiation  toughness  or  the  resistance 
curve  shape. 

O’Dowd  aid  Shih  (1991)  have  proposed  that  the  difference  in  crack  tip  stress  fields 
between  the  small  scale  yielding  (SSY)  finite  element  solution  and  the  finite  body  solution  can 
be  quantified  in  terms  of  a  field  quantity  that  they  have  call  Q.  The  Q  quantity  is  a  function  of 
J,  the  crack  shape  and  size,  the  structural  geometry,  mode  of  loading  and  on  the  level  of 
deformation  and  can  only  be  calculated  from  a  high  resolution  elastic-plastic  computational 
analysis.  O’Dowd  and  Shih  propose  that  a  J-Q  fracture  locus  can  be  developed  experimentally 
for  a  particular  material,  with  higher  Q’s  meaning  higher  constraint  and  resulting  in  lower  J 
measures  at  cleavage  or  at  ductile  crack  initiation,  or  at  whatever  critical  measurement  point  is 
to  be  used.  This  procedure  avoids  the  need  to  start  with  a  fracture  criterion,  which  is  a  serious 
difficulty  if  one  wishes  to  apply  the  Dodds  and  Anderson  approach  to  the  case  of  ductile  crack 
extension. 

A  similar,  simpler,  but  more  controversial  approach  has  been  suggested  by  Bet6gon  and 
Hancock  (1991),  who  use  the  non-singular  term  of  the  elastic,  Williams  (1957),  crack  singularity 
solution,  called  the  T-Stress",  as  a  measure  of  elastic-plastic  crack  tip  constraint.  This  quantity 
is  only  dependent  on  the  initial  geometry  and  loading  and  hence  is  relatively  easy  to  calculate, 
and  is  available  in  the  literature  for  many  geometries. 

The  objective  of  this  work  is  to  develop  some  upper  shelf,  elastic-plastic  experimental 
results  to  attempt  to  investigate  the  applicability  of  the  Q  and  T  stress  parameters  to  the 
correlation  of  upper  shelf  initiation  touglmess  and  J  resistance  curves.  The  fu^t  objective  was 
to  obtain  upper  shelf  J  resistance  curves,  1,^,  and  tearing  resistance,  T„„„  results  for  a  range  of 
applied  constraint.  The  J-Q  and  J-T  stress  loci  were  developed  and  compared  with  the 
expectations  of  the  O’Dowd  and  Shih  and  the  Betegon  and  Hancock  analyses.  Constraint  was 
varied  by  changing  the  crack  length  and  also  by  changing  the  mode  of  loading  from  bending  to 
predominantly  tensile. 

The  principal  conclusions  of  this  work  are  that  Ji^  does  not  appear  to  be  dependent  on  T 
stress  or  Q  while  the  material  tearing  resistance  is  dependent  on  T  stress  and  Q,  with  the  tearing 
modulus  increasing  as  constraint  decreases. 
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OB.TECTrVES 


Constraint  has  been  an  important  consideration  in  fracture  mechanics  from  the  earliest 
work  that  was  done  to  develop  the  1974  version  of  the  ASTM  Standard  E399.  Stringent 
thickness  and  size  requirements  were  placed  on  the  test  samples,  in  terms  of  the  measured  plastic 
zone  size.  These  requirements  often  invalidated  the  results  in  practical  application  sizes.  The 
restrictions  have  been  retained  in  later  revisions,  often  requiring  the  engineer  to  test  in  his 
application  thickness  and  to  work  generally  with  “invalid"  data,  using  engineering  judgement  and 
experience  to  avoid  catastrophic  failure  of  his  structure.  Similar  size  requirements  have  been 
incorporated  into  mote  recent  elastic-plastic  fracture  test  standards  like  ASTM  E813  and  ASTM 
El  152,  even  though  experimental  data  has  not  shown  a  clear  and  direct  relationship  between 
specimen  size  and  the  resulting  initiation  toughness  or  the  resistance  curve  shape. 

Recent  computational  work  by  Anderson  and  Dodds  (1991)  and  Dodds,  Anderson  and 
Kirk  (1991)  has  sought  to  quantify  constraint  by  comparing  the  stress  field  near  a  crack  in  an 
elastic-plastic  body  to  the  small-scale  yielding  (SSY)  solution  that  corresponds  to  that  geometry 
and  material  combination,  as  also  determined  by  numerical  computation.  This  SSY  solution  has 
been  shown  to  be  basically  different  than  the  Hutchinson  (1968),  Rice  and  Rosengren  (1968) 
(HRR)  solution,  apparendy  due  to  the  crack  tip  blunting  which  is  invariably  present  near  the 
crack  dp.  The  Dodds,  et  al  work  has  shown  that  constraint  depends  on  specimen  or  structure 
thickness,  on  in-plane  dimensions,  and  on  the  mode  of  loading.  They  also  showed  that  the 
principal  stress  fields  in  two-dimensional  cases  are  self-similar  and  a  comparison  of  similar 
stressed  volumes  can  be  made  which  allows  a  comparison  of  constraint  from  one  situation  to 
another,  i.e.  a  short-cracked  bend  geometry  can  be  chosen  to  match  the  constraint  of  a  much 
harder  to  test  shon  cracked  tensile  geometry.  In  addition,  they  developed  a  scaling  model  to 
predict  the  effect  of  changes  in  constraint  on  the  cleavage  h'acture  toughness. 

The  Dodds,  et  al.  scaling  model  has  been  applied  predominantly  to  cleavage  fracture, 
since  as  presently  formulated,  it  assumes  a  fracture  criteria  similar  to  that  proposed  by  Ritchie, 
Knott,  and  Rice  (1973)  (fracture  occurs  when  the  maximum  principal  stress  reaches  a  critical 
value.  Of,  at  a  critical  distance,  r*,  from  the  crack  tip).  By  scaling  the  remote  J  integral  to  achieve 
similar  stress  conditions  near  the  crack  tip  a  correction  of  the  J  integral  values  at  cleavage 
initiation  can  be  made  which  accounts,  at  least  approximately,  for  the  effects  of  constraint  on 
fracture.  The  applicability  of  this  technique  in  the  lower  shelf  ductile-brittle  transition  regime 
of  ferritic  steels  has  been  demonstrated  by  Sorem,  Dodds  and  Rolfe  (1991)  and  by  Kirk, 
Koppenhoefer  and  Shih  (1993). 

O’ Dowd  and  Shih  (1991)  have  proposed  that  the  difference  in  crack  tip  stress  fields 
between  the  SSY  solution  and  the  finite  body  solution  can  be  quantified  in  terms  of  a  field 
quantity  that  they  have  call  Q.  The  Q  quantity  is  a  function  of  J,  the  crack  shape  and  size,  the 
structural  geometry,  mode  of  loading  and  on  the  level  of  deformation  and  can  only  be  calculated 
from  a  high  resolution  elastic-plastic  computational  analysis.  O’Dowd  and  Shih  propose  that  a 
J-Q  fracture  locus  can  be  developed  experimentally  for  a  particular  material,  with  higher  Q’s 
meaning  higher  constraint  and  resulting  in  lower  J  measures  at  cleavage  or  at  ductile  crack 
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initiation,  or  at  whatever  critical  measurement  point  is  to  be  used.  This  procedure  avoids  the 
need  to  start  with  a  fracture  criterion,  which  is  a  serious  difficulty  if  one  wishes  to  apply  the 
Dodds  and  Anderson  approach  to  the  case  of  ductile  crack  extension.  Examples  of  a  J-Q  fracture 
locus  for  the  case  of  cleavage  have  been  developed  by  Kirk,  Koppenhoeffer  and  Shih  (1993)  and 
also  by  Sumpter  and  Forbes  (1993),  and  these  results  seem  to  support  the  J-Q  concept. 

A  similar,  simpler,  but  more  controversial  approach  has  been  suggested  by  Betfgon  and 
Hancock  (1991),  who  use  the  Erst  non-singular  term  of  the  elastic  crack-tip  stress  field  solution 
(Williams  1957)  as  a  measure  of  elastic-plastic  crack  tip  constraint.  This  quantity  has  been  called 
the  'T-stress",  a  notation  used  by  Larsson  and  Carlsson  (1973)  who  studied  the  effect  of  this  term 
on  the  crack  tip  stress  Eeld.  This  quantity  is  only  dependent  on  the  initial  geometry  and  loading 
and  hence  is  relatively  easy  to  calculate,  and  is  available  in  the  literature  for  many  geometries. 
The  application  of  this  elastic  quantity  to  elastic-plastic  fracture  is  highly  controversial,  but 
published  results  have  shown  good  correlation,  and  the  ease  of  use  is  a  strong  selling  point  for 
this  parameter. 

The  objective  of  this  work  is  to  develop  some  upper  shelf,  elastic-plastic  experimental 
results  to  investigate  the  applicability  of  the  Q  and  T  stress  parameters  to  the  correlation  of  upper 
shelf  initiation  toughness  and  J  resistance  curves.  The  Erst  objective  was  to  obtain  upper  shelf 
J  resistance  curves,  J,^,  and  tearing  resistance,  T„„(,  results  for  a  range  of  applied  constraint.  The 
J-Q  and  J-T  stress  loci  were  developed  and  compared  with  the  expectations  of  the  O’Dowd  and 
Shih  and  the  Betegon  and  Hancock  analyses.  Constraint  was  varied  by  changing  the  crack  length 
and  also  by  changing  the  mode  of  loading  from  bending  to  pr^ominantly  tensile.  Test 
techniques  and  analysis  have  been  developed  as  needed  for  the  low  constraint  fracture  test 
geometries.  Two  materials  have  been  used  in  this  study,  an  HY-100  high  strength  structural  steel 
and  an  A533B  pressure  vessel  steel.  Some  of  the  results  for  the  HY-100  steel  have  been  reported 
previously  in  NUREG/CR-5879  (Joyce,  et  al.  1992),  they  are  repeated  again  here  to  demonstrate 
the  consistency  of  the  results  found  for  the  two  materials. 
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1.0  EXPERIMENTAL  DETAILS 


1.1  Material  Description 

Two  structural  steels  were  tested  in  this  study.  The  Hrst  material  was  HY-100,  a  high 
strength  structural  steel  with  tensile  mechanical  properties  and  chemistry  as  shown  in  Table  1. 
This  material  was  from  a  6.35  cm  thick  plate  and  all  specimens  were  oriented  so  that  the  crack 
plane  was  in  the  T-L  direction  as  designated  by  ASTM  E399.  The  second  material  was  an 
ASTM  A533  Grade  B  pressure  vessel  steel  with  the  tensile  mechanical  properties  and  chemistry 
also  shown  in  Table  1.  The  plate  for  this  material  was  originally  21.5  cm  thick,  but  for  these 
tests  all  samples  were  cut  from  the  center  15  cm.  All  specimens  of  this  material  were  oriented 
in  the  L-T  orientation. 

1.2  Specimen  Details 

Four  distinct  test  geometries  were  studied  in  this  work:  the  standard  IT  compact  specimen 
C(T),  the  single  edge -notched  bend  specimen  SE(B),  the  single  edge-notched  tensile  specimen, 
SE(T)  that  was  pin-loaded  and  the  double  edge-notched  tensile,  DE(T),  specimen,  also  pin-loaded. 
The  SE(B)  specimens  were  tested  in  a  standard  deep  notched  configuration  with  crack  length 
specimen  width  ratio  a/W  -  0.6,  and  also  in  a  shallow  notched  configuration  with  a/W  -  0.15. 
TTie  SE(T)  and  DE(T)  conHgurations  were  tested  with  a/W  ratios  from  0.35  to  0.7.  Schematic 
drawings  of  the  SE(T)  and  DE(T)  geometries  are  shown  in  Figure  1 .  All  specimens  were  25  mm 
thick  and  side  grooved  to  a  total  thickness  reduction  of  20%.  All  specimens  were  side  grooved 
after  precracking. 

Some  non-side  grooved  specimens  of  the  HY-100  steel  were  tested  and  reported  in  Joyce, 
et  al.  (1992),  however  sdl  sfiecimen  results  presented  here  are  for  specimens  which  have  been  side 
grooved  to  a  total  reduction  of  20%.  Matrices  of  the  test  specimens  are  presented  in  Table  2  and 
3  for  the  two  materials  tested  and  repiorted  here. 

Tests  of  the  HY-100  steel  were  done  at  ambient  temperature  (25°C),  wnile  the  A533B  was 
tested  at  approximately  100°C  to  assure  that  the  fracture  mode  for  all  tests  were  fully  ductile 
throughout. 

1.3  Specimen  Precracking 

The  bend  specimens  and  the  tensile  bars  were  precracked  in  bending  using  a  three  ]X)int 
bend  apparatus.  The  short  cracked  HY-100  bend  specimens  were  precracked  starting  from  a  wide 
specimen,  with  W  -  70  mm,  and  precracked  until  the  crack  was  about  27  mm  long.  The 
specimens  were  subsequently  machined  to  remove  the  material  at  the  crack  flanks,  until  a  final 
configuration  was  obtained  with  a  crack  length  of  about  7  mm  in  a  remaining  ligament  of  50 
mm.  The  precrack  fronts  obtained  in  this  fashion  for  this  material  were  found  to  be  straight  and 
accurate  in  all  cases,  but  the  method  was  expensive  and  arduous.  The  A533B  short  crack  bend 
specimens  were  precracked  by  starting  with  a  bar  of  width  50  mm  with  a  short  machined  notch 
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Figure  1  Schematic  drawings  of  the  SE(T)  and  DE(T)  geometries  tested  in  this 
investigation. 
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S  mm  in  depth.  Then  a  single  large  reverse  bending  load  was  applied  that  was  calculated  to 
introduce  a  compressive  plastic  zone  of  0.S  mm  in  extent.  The  specimen  was  then  fatigue 
precracked  as  usual.  The  crack  length  was  monitored  using  a  computer  controlled  servo- 
hydraulic  test  machine,  and  the  short  crack  compliance  equation  obtained  by  Joyce  (1992),  until 
the  final  precrack  length  of  approximately  7  mm  was  obtained.  These  cracks  were  found  to  be 
straight,  and  the  single  load  reversal  method  was  much  simpler  than  the  double  machining 
process  used  on  the  HY-100  short  cracked  bend  specimens.  This  technique  has  become  standard 
practice,  even  on  large,  shallow  cracked  bend  bars. 

The  SE(T)  specimens  were  precracked  in  three  point  bending,  starting  with  machined 
notches  with  a/W  -  0.15,  and  grown  to  a/W  values  of  between  0.35  and  0.65  for  testing.  The 
DE(D  specimens  were  precracked  in  four  point  bending,  with  two  rollers  placed  tightly  across 
the  compression  side  notch.  The  initial  notch  depths  were  kept  at  a/W  -  0.15  in  the  DE(T) 
specimens  even  though  deep  cracks  were  desired  so  that  the  ligament  during  precracking  was  a 
large  as  possible.  A  tightly  fitting  wedge  was  pushed  into  the  compression  side  notch,  and  this 
allowed  using  standard  bend  compliance  equations  for  monitoring  the  crack  length  of  the  DE(T) 
specimens  during  precracking.  The  specimens  were  reversed  several  times  to  obtain  even 
precracks  on  both  sides.  Matching  the  bend  compliances  in  the  two  directions  seemed  to 
accurately  match  the  lengths  of  the  two  cracks.  As  a  final  check  the  specimen  was  loaded  in 
tension  and  the  outputs  of  two  clip  gages  mounted  across  the  two  cracks  were  compared.  It  was 
generally  found  that  the  bending  compliance  matched  cracks  were  of  equal  length  and  any 
remaining  difference  in  COD  output  was  usually  do  to  misalignment  of  the  test  machine  load 
train. 

1 .4  Test  Technique 

All  tests  were  conducted  using  a  single  specimen,  computer  interactive,  unloading 
compliance  test  procedure  which  allowed  monitoring  the  specimen  crack  length  and  the  applied 
J  integral  during  the  course  of  the  test.  Equations  were  presented  in  Joyce,  et  al.  (1992)  for  the 
required  elastic  and  plastic  components  of  the  J  integral  and  for  estimating  the  crack  length  or 
lengths  from  the  experimentally  measured  crack  opening  compliance.  A  simple  rotation 
correction  procedure  greatly  improves  the  crack  length  estimation  accuracy  of  the  SE(T) 
specimen.  The  rotation  correction  that  is  used  here  is  similar  to  that  used  in  ASTM  El  152  for 
the  C(T)  specimen,  a  correction  that  is,  of  course,  used  for  the  C(T)  specimens  tested  in  this 
work.  The  SE(T)  specimen  rotation  correction  is  developed  in  the  next  section.  The  results 
presented  here  for  the  HY-100  alloy  SE(T)  specimens  are  altered  considerably  from  that  reported 
previously  (Joyce,  et  al.  1993)  because  the  SE(T)  specimen  rotation  correction  is  used  here. 

In  all  cases,  crack  growth  corrected  J  equations  are  used  here,  as  developed  in  Joyce,  et 
al.  (1992),  which  are  similar  to  those  required  for  the  deep  SE(B)  specimens  and  the  C(T) 
specimens  by  ASTM  El  152.  All  data  was  stored  on  magnetic  media  for  subsequent  re-analysis 
as  needed. 

The  SE(B)  specimens  were  tested  with  standard  bend  fixtures  which  were  made  much 
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Table  2 


List  of  HY-100  specimens  tested  in  this  investigation. 


Specimen 

ID 

Type 

a/W 

FYO  1 

SE(B) 

0.66 

FY0  3 

SE(B) 

0.66 

FYO  21 

SE(B) 

0.14 

FYO  26 

0.13 

FYO  27 

SE(B) 

0.14 

FYO  150 

SE(B) 

0.61 

FYO  151 

SE(B) 

0.61 

FYO  158 

SE(B) 

0.60 

FYO  159 

SE(B) 

0.62 

FYO  160 

SE(B) 

0.11 

FYO  161 

SE(B) 

0.11 

FYO  2SB 

SECT) 

0.40 

FY0  3SB 

SECT) 

0.47 

FY0  4SA 

SECT) 

0.65 

FYO  lOSA 

SE(T) 

0.35 

FYO  USB 

DE(T) 

0.68 

FYO  12SA 

DE(T) 

0.61 

B 

Bn 

W 

(mm) 

(mm) 

(mm) 

50. 

40. 

50. 

50. 

40. 

50. 

50. 

40. 

50. 

25. 

20. 

50. 

25. 

20. 

50. 

25. 

20. 

50. 

25. 

20. 

50. 

12.5 

10. 

50. 

12.5 

10. 

50. 

12.5 

10. 

50. 

12.5 

10. 

50. 

25. 

20. 

64. 

25. 

20. 

64. 

25. 

20. 

64. 

25. 

20. 

64. 

25. 

20. 

32. 

25. 

20. 

32. 

Side 

Groove 

(Y/N) 


Table  3  List  of  AS33B  specimens  tested  in  this  investigation. 


Specimen  ID 

Type 

a/W 

Side 

Groove? 

(Y/N) 

B 

(mm) 

Bn 

(mm) 

W 

(mm) 

cr3 

C(D 

0.6 

Y 

25. 

20. 

50. 

CT9 

C(T) 

0.6 

Y 

25. 

20. 

50. 

CTIO 

C(T) 

0.6 

Y 

25. 

20. 

50. 

DBl 

SE(B) 

0.62 

Y 

25. 

20. 

50. 

DB2 

SE(B) 

0.62 

Y 

25. 

20. 

50. 

DB3 

SE(B) 

0.62 

Y 

25. 

20. 

50. 

SBl 

SE(B) 

0.15 

Y 

25. 

20. 

50. 

SB2 

SE(B) 

0.15 

Y 

25. 

20. 

50. 

SB3 

SE(B) 

0.15 

Y 

25. 

20. 

50. 

SENl 

SECT) 

0.41 

Y 

25. 

20. 

63.5 

SEN2 

SECT) 

0.38 

Y 

25. 

20. 

63.5 

SEN4 

SECT) 

0.66 

Y 

25. 

20. 

63.5 

SEN9 

SECT) 

0.41 

Y 

25. 

20. 

63.5 

SENIO 

SE(T) 

0.62 

Y 

25. 

20. 

63.5 

SE5D 

DECT) 

0.7 

Y 

25. 

20. 

31.8 

SE6D 

DE(T) 

0.7 

Y 

25. 

20. 

31.8 

SE7D 

DE(D 

0.68 

Y 

25. 

20. 

31.8 
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sturdier  than  usual  to  accommodate  the  higher  loads  typical  of  short  crack  specimens.  The 
procedures  used  for  the  short  crack  tests  were  presented  by  Joyce  (1992)  which  showed  clearly 
that  unloading  compliance  was  a  viable  test  technique  for  SE(B)  specimens  with  aAV  as  short 
as  0.1.  The  compliance  equations  are  much  less  sensitive  in  the  short  crack  region,  but  the  load 
applied  increases  with  (W-a)^  so  that  for  the  short  crack  specimens  the  unloadings  become  much 
larger,  and  the  crack  opening  displacement  (COD)  continues  to  be  adequately  large  and  can  be 
measured  with  a  high  resolution  digital  voltmeter.  The  combination  of  high  loads  and  less 
sensitive  unloading  compliance  makes  these  test  more  difficult,  but  if  care  is  taken,  excellent 
results  can  be  obtained.  A  flex  bar  was  used  to  measure  the  load  line  displacements  for  all  SE(B) 
specimens  as  described  previously  (Hacked  and  Joyce,  1986)  since  significant  indentations  did 
occur  at  the  rollers  for  these  specimens. 

The  SE(T)  and  DE(T)  specimens  were  loaded  with  oversized  tension  clevises  similar  to 
what  is  used  for  standard  C(T)  specimens.  The  HY-100  tests  were  done  with  clevises  that  had 
round  holes  while  the  A533B  tests  were  done  with  clevises  that  had  holes  with  loading  flats  to 
allow  for  the  loading  pin  to  roll  as  the  specimen  rotated  during  test.  It  was  found  that  allowing 
the  specimen  half  to  rotate  freely,  and  correcting  for  the  rotation  effects  was  the  preferred  method 
for  the  SE(T)  tests.  For  DE(T)  specimens,  however,  it  was  found  that  round  holes  were 
preferred,  providing  an  initial  alignment  that  was  essential  to  accurately  test  the  DE(T)  specimen 
both  at  the  start  of  test  and  as  crack  growth  proceeded.  For  the  SE(T)  specimens  a  standard  clip 
gage  was  installed  to  measure  the  crack  mouth  opening  displacement,  which  was  used  for  crack 
length  estimation,  and  an  LVDT  gage  was  installed  on  the  initial  specimen  load  line  to  measure 
the  load  line  extension  of  the  specimen.  For  the  DE(T)  specimen,  two  COD  gages  were  used 
as  well  as  an  LVDT  gage  on  the  specimen  centerline.  In  general  the  average  COD  displacement 
was  used  to  estimate  the  average  crack  length  for  the  DE(T)  specimen.  Both  COD  gage  readings 
were  recorded  in  the  data  file  and  can  be  plotted  separately,  if  desired. 
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2.0  ANALYSIS 
2.1  J  Intearal  Analysis 

The  J  integral  is  calculated  here  by  summing  the  elastic  and  plastic  components,  with  the 
components  calculated  separately.  The  elastic  J  component,  Jd,  is  calculated  ^m; 


(1) 


where  K  is  the  elastic  stress  intensity  factor  for  the  specimen,  E'-E/(l-v^),  and  E  and  v  are  the 
elastic  modulus  and  Poisson’s  ratio,  respectively.  The  plastic  J  component,  Jpj,  is  calculated  using 
the  ASTM  Standard  El  152  Jp,  equation: 


•Vo  *  "V*-!)  * 


Hi 


Uo  ~  V<-i) 


s 


1  - 


Y,  (fly  -  fl(y-i)) 


(2) 


with: 

Apu  -  area  under  the  load  versus  plastic  load  line  displacement  curve  to  increment  i, 
B,4  -  net  specimen  thickness  at  the  side  groove  roots, 

Tii  -  the  plastic  T]  factor  at  crack  length  a^ 
bj  -  the  incremental  remaining  ligament 
W  -  the  specimen  width  and 


Yy  = 


(3) 


evaluated  at  the  crack  length  aj,  and: 


Ti/  *  dr\J<KaJW) 


(4) 


Formulas  for  the  K’s,  T|’s,  and  y’s  used  for  the  SE(B),  SE(T),  and  DE(T)  specimens  are  presented 
in  the  next  subsections. 

2.2  SE(B)  Analysis 

Previous  work  by  Joyce  (1992)  has  shown  that  unloading  compliance  can  be  used  to 
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evaluate  J-R  curves  for  short  crack  bend  specimens.  As  the  crack  becomes  very  short  the 
compliance  equation  becomes  less  sensitive  to  crack  length  but  the  specimen  limit  load  also 
increases,  which  increases  the  length  of  the  allowed  elastic  unloading,  and  the  total  effective 
crack  length  measurement  resolution  is  only  slightly  degraded.  Results  obtained  by  Joyce  (1992) 
speared  to  be  fully  adequate  for  J,c  and  J-R  curve  testing  for  a/W  ratios  as  small  as  0.15.  In 
this  work  similar  success  was  found  for  a/W  ratios  as  small  as  0.1.  To  test  in  this  a/W  range  a 
new  equation  to  estimate  crack  length  from  the  specimen  COD  compliance  is  needed  since  the 
equation  available  in  ASTM  E813  and  El  152  does  not  apply  for  a/W  ratios  below  0.4.  An 
equation  for  bend  specimen  compliance  as  a  function  of  a/W  that  is  good  for  all  a/W  is  available 
in  The  Stress  Analysis  of  Cracks  Handbook  (The  Handbook)  (Tada,  Paris  and  Irwin  1985).  This 
equation  is: 


P 


_  24(alW) 


BWE' 


I  S/4 


0.76  -  2.28 

f  "  ]  +  3.87f  "  f  -  2.04f  " 

[w]  [wj  (l-a/W)\ 

(5) 


where:  5  -  crack  mouth  opening  displacement  at  the  specimen  edge 

P  -  load 

B  -  specimen  thickness 
S  -  specimen  span 

For  the  short  crack  range  a  reverse  fit  to  calculate  the  crack  size,  a/W,  from  the  measured 
compliance  was  obtained  by  Joyce  (1992)  using  a  standard  fifth  order  polynomial  by  restricting 
the  a/W  range  to  between  0.05  and  0.45.  The  fit  is  accurate  to  within  0.06%  which  is  acceptable 
for  the  unloading  compliance  method.  The  Joyce  relationship  is: 

—  =  1.01878  -  4.5367tt  +  9.0101  -  27.333u^  +  74.4«^  -  71.489ii® 

W 


u 


1 


I  S/4 


+  1 


(6) 


where  C  -  unloading  compliance,  S/P,  and  has  been  used  for  the  short  crack  SE(B)  specimens 
presented  below.  The  standard  equation  of  ASTM  El  152  was  used  for  the  deep  cracked  bend 
specimens  analyzed  below. 

For  the  deep  cracked  SE(B)  specimens  the  q  and  y  factors  of  ASTM  El  152,  (q  -  2.0  and 
y  -  1.0),  are  used  in  Eq.  2  to  evaluate  J.  For  the  short  crack  specimens,  however,  these 
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coefficients  must  be  changed  to  accurately  evaluate  J.  This  problem  has  been  investigated  by 
Haigh  and  Richards  (1974),  Sumpter  (1987),  and  by  Joyce  (1992).  A  comparison  of  various 
estimates  of  T]  is  shown  in  Figure  2  which  includes  results  of  the  above  authors  and  results 
derived  by  Joyce  (1992)  from  Elastic-Plastic  Fracture  Analysis  (EPRI  Handbook)  (Kumar, 
German  and  Shih  1981).  The  ABAQUS  results  were  obtained  by  Joyce  (1992)  using  a  2D  elastic 
plastic  finite  element  analysis  incorporating  incremental  plasticity.  In  the  work  that  follows  the 
polynomial  function  for  iii  developed  by  Sumpter  (1987)  is  us^  for  all  short  cracked  SE(B) 
specimens  with  a/W  <  0.282.  This  polynomial  expression  is: 


Ti  =  0.32  +  ilifllW)  -  49.S{alW)^  +  99.8(fl/»0’ 

This  equation  gives  ■q  <  2.0  for  aAV  <  0.282.  Sumpter  switches  to  q  -  2,0  when  the  crack 
length  exceeds  a/W  -  0.282.  In  this  work  the  short  crack  specimens  were  started  and  completed 
with  a/W  <  0.282.  Equation  (7)  was  used  because  it  was  complete  in  the  range  of  interest  and 
because  of  the  correspondence  with  the  EPRI  Handbook  results. 

The  7  factor  is  calculated  from  q  using  Eq.  3.  For  the  short  crack  specimens  y  was 
obtained  by  differentiating  Eq.  (7)  to  give: 


-12.22 +  106.7| 

'a' 

-236.6 

[-1 

^-924.6| 

'a' 

k 

j^+4845.4| 

a' 

k 

*-9880 

’+9960 

f-t 

0.32  +  12 

(-1 

-49.5j 

a' 

^  +99.8 

\w] 

|3 

(8) 


Recent  work  by  Kirk  and  Dodds  (1993)  has  suggested  that  J  should  be  calculated  using 
an  qcoD  b^ised  on  the  crack  mouth  opening  displacement  rather  than  on  the  load  line 
displacement.  Kirk  and  Dodds  show  with  computational  calculations  that  q(^Q  is  much  less 
dependent  on  strain  hardening  than  is  the  standard  load  line  q.  The  equation  used  then  to 
calculate  J  using  q^oo  is: 


^COD^CODpKO 


b,B, 


N 


(9) 


where:  ^codpki)  ~  under  the  load  versus  plastic  crack  opening  displacement  curve 

to  increment  i. 

This  equation  does  not  include  the  crack  growth  correction  term  used  in  all  other  equations  in 
this  work,  and  this  causes  it  to  have  serious  problems  since  the  crack  growth  correction  is 
essential  for  the  development  of  a  correct  J-R  curve  beyond  1  mm  of  crack  growth. 
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Predicted  plastic  “n  factors  for  SE(B)  specimens. 


For  unloading  compliance  testing  of  SE(T)  specimens  equations  are  required  for  K,  t),  and 
y,  as  functions  of  a/W.  The  equations  used  in  this  work  are  presented  in  the  following  sections. 


2.3.1  SECn  K  Expression 


Since  the  SE(T)  specimens  tested  here  had  pin  loading,  the  K  expressions  for  fixed  end 
loading  in  The  Handbook  were  checked  with  ABAQUS  finite  element  analysis.  A  total  of  14 
different  SE(T)  finite  element  grids  were  developed  with  0.12  <  a/W  <  0.80.  These  grids  were 
used  to  develop  both  the  elastic  stress  intensity  factor  K  and  the  plastic  t|  factor  as  described 
below.  The  stress  intensity  factor  relationship  was  assumed  to  have  the  form: 


(10) 


and  F(a/W)  was  fit  with  a  polynomial  to  give: 

F(a/W)  =  -0.0917  +  22.392(fllW)  -  141.96(q/R0^  +  449.72(a/»0’ 

-  M5.59{alWf  +  363.52 (n/HO’ 


(11) 


This  equation  fit  the  ABAQUS  results  within  ±2%  over  the  a/W  range  from  0.12  to  0.80. 
A  comparison  is  presented  in  Figure  3  with  the  ABAQUS  results  and  a  standard  form  taken  from 
The  Handbook.  Qearly  The  Handbook  equation,  the  polynomial  fit,  and  the  ABAQUS  results 
agree  very  well  in  the  range  of  0.12  <  a/W  <  0.80.  In  the  experimental  work  presented  below 
the  polynomial  form  for  F(a/W)  presented  in  Eq.  (10)  has  been  used  for  all  SE(T)  specimens. 


2.3.2  SEm  n  Factor 

In  Joyce,  et  al.  (1992)  several  methods  were  used  to  estimate  the  “n  foi*  Oie  SE(T) 
specimen  including  elastic-plastic  finite  elements  analysis,  the  EPRI  Handbook,  and  published 
results  by  Wu,  et  al.  (1990)  and  Sharobeam,  et  al.(1991).  A  comparison  of  all  of  these  results 
from  Joyce,  et  al.  (1992)  is  shown  in  Figure  4. 


A  ABACUS  -  KYI 00 
0  SX-WU  (Modified  EFRI) 
o  SX-WU  Limit  Load 
0  EPRI  Plane  Strain  n  • 
o  EPRI  Plane  Stress  n  « 
I - 1  Sharobeam  and  Landes 


5  0.6  0.7  0.8  0.9 


SECT)  geometry 


In  the  experimental  work  that  follows  the  dashed  bi-linear  relationship  shown  in  Figure  4 
was  used  to  evaluate  tIj  at  each  crack  length  a^.  This  form  also  allowed  calculating  y.  front  Eq. 
(3)  which  is  necessary  to  calculate  Jp,  using  Eq.  (2).  The  equations  used  to  evaluate  Tij  and  y  for 
the  SE(T)  specimen  are  thus: 


Ti,  =  5.71  (flJW) 

0  a JW  ^  0.417 

(12) 

Ti^  =  2.38 

0.417  <  aJW  i  1.0 

(13) 

(5  71) 

-l-ibJW)  ^ 

0  <  ajW  i  0.417 

(14) 

Y,  =  1.38 

0.417  <  aJW  i  1.0 

(15) 

2.3.3  SETT)  Crack  Length  Estimation 


Since  the  SE(T)  specimen  is  of  a  rather  short  length  and  has  the  load  applied  through  the 
centered  pin  holes,  the  compliance  equations  in  standard  fracture  mechanics  handbooks  like  The 
Handbook  are  not  necessarily  applicable.  The  standard  forms  available  assume  uniform  stresses 
at  the  loading  edges  and  the  SE(T)  configuration  used  here  was  not  thought  to  be  long  enough 
to  allow  the  direct  use  of  equations  based  on  the  uniform  stress  assumption.  A  finite  element 
analysis  was  used  in  Joyce,  et  al.  (1992)  to  develop  a  polynomial  equation  giving  the  crack 
length  as  a  function  of  the  COD  compliance  for  the  SECT ,  specimen  geometry  used  here.  This 
equation  has  the  form: 

alW  =  1 .012525  -  2.95323(«  V  +  6.68(h  f  -  17.7954(tt  f  (16) 

+ 25.35 17(tty-12.9747(uy 


with: 


u'  = 


1  + 


N 


E'Bh 


(17) 
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For  side  grooved  specimens  the  thickness  B  is  replaced  by  where: 


B 


(18) 


where  B^  is  the  net  specimen  thickness  at  the  side  groove  roots.  This  effective  thickness 
formulation  is  consistent  with  ASTM  E813  and  El  152. 
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2.3.4  SEfD  Rotation  Correction 


A  rotation  correction  can  be  developed  for  the  SE(T)  specimen  using  the  notation  of 
Figure  5.  Two  separate  corrections  are  needed,  one  to  correct  the  knife  edge  COD  displacement 
for  the  effect  of  rotation,  the  second  to  correct  the  load  for  the  effect  of  rotation  (Loss  1977; 
Merkle,  198x). 

The  objective  of  the  first  correction  is  to  correct  the  measured  displacement,  Ad^ 
increment  during  an  unloading  to  obtain  the  corrected  displacement  increment,  Ad^^,  as  shown 
in  Figure  5a.  Using  the  geometry  of  Figure  5a,  the  rotation  angle,  0,  is: 


0  = 


id^  -  D) 


]  -  tan'* 


_D 


(19) 


Using  similar  triangles,  it  can  be  shown  that : 

^ 

Ad^  R-y 


(20) 


where  y  is  the  horizontal  shift  of  the  displacement  measurement  point  due  to  rotation  of  the 
specimen.  The  quantity  (R-y)  is  obtained  from  the  relationship 


cos(0  +  P) 


R-y 
*  R} 


(21) 


Using  the  identity  for  cos(x+y) 

cos(6  +  P)  =  cosScosp  -  sinOsinp 


(22) 


and  noting  that 


sinP  = 


D 


cosp  =  ^ 


(23) 


eq.  (21)  can  be  solved  for  (R-y).  Substituting  the  resulting  expression  for  (R-y)  into  (20)  and 
rearranging  terms  leads  to  the  following  expression 


(24) 
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Figure  5  Geometric  relationship  of  variables  used  in  rotation  correction  development. 
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This  equation  was  used  to  correct  the  measured  load  line  compliance  Ci^Dm  to  obtain  the 
corrected  load  line  compliance,  Ci,rv,  before  calculating  the  elastic  and  plastic  area  components 
used  to  calculated  the  elastic  and  plastic  J  components. 

To  use  this  aiudysis  it  is  necessary  to  assume  a  position  for  the  center  of  rotation  for  the 
SE(T)  specimen.  In  the  standard  C(T)  analysis  of  ASTM  El  152  the  center  of  rotation  is  assumed 
to  be  at  the  center  of  the  remaining  ligament,  and  that  assumption  was  also  used  here,  i.e.  - 
(a  +  Wy2  and  Rl  -  Rg  •  W/2. 
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For  unloading  compliance  testing  of  the  DE(T)  specimens,  equations  are  required  for 
calculation  of  K,  ii,  and  y  as  functions  of  a/W  and  for  a/W  as  a  function  of  the  COD  compliance, 
fi/P.  The  equations  used  for  this  projea  are  presented  in  the  following  sections. 

2.4.1  DETD  K  Expression 

The  K  equation  for  the  DE(T)  specimen  with  a  deeply  cracked  geometry  can  be  taken 
directly  from  The  Handbook.  The  equation  used  has  the  form; 


P 

2WB 


Fia/m 


(30) 


with: 


F(a/W) 


1.122 -0.561  f— 1+0.205 

_ [w 


[— f  +0.471  f—f  -O.lof— f 

[wj  [  lyj  i  lyj 


(31) 


yjl-alW 


This  equation  should  be  accurate  to  ±0.5%  for  any  a/W,  but  is  limited  to  a/W  >  0.6  by  the  pin 
hole  loading.  Three  finite  element  computations  by  Joyce  et  al.(1992)  showed  that  this  equation 
is  accurate  within  ±2%  for  deeply  cracked  DE(T)  specimens  including  the  center  pin  holes. 

2.4.2  DE(T)  n  and  v  Factors 

The  T)  factor  for  the  DE(T)  specimen  geometry  was  obtained  from  both  elastic-plastic 
finite  element  analysis  using  ABAQUS  and  from  the  EPRI  Handbook.  The  T|  factor  used  here 
is  taken  to  relate  the  J  integral  at  each  crack  to  the  total  plastic  work  applied  to  the  specimen, 
i.e.: 


V 


(32) 


where: 

Ap,  -  plastic  area  under  the  specimen  load  versus  plastic  load  line  displacement  plot 
b  -  specimen  half  remaining  ligament 
B  -  specimen  thickness 
1)  -  plastic  Ti  factor. 


Analytical  work  by  Wu,  et  al.  (1990)  based  on  limit  load  theory,  shows  that  should  be 


nearly  constant  for  the  DE(T)  specimen  over  the  a/W  range  of  interest  here.  A  value  of 
q)proximately  0.27,  instead  of  the  usual  2.0,  is  also  predicted  in  Wu,  et  al.  (1990)  with  only  a 
very  slight  dependence  on  strain  hardening.  These  predictions  are  confirmed  here  by  both  the 
finite  element  analysis  and  the  EPRI  analysis. 

Deeply  notched  DE(T)  specimens  with  three  different  crack  sizes  were  analyzed  using 
finite  elements  and  the  results  shown  in  Figure  6,  compared  to  the  results  of  Wu,  et  al.  (1990) 
Other  results  shown  on  Figure  6  include  calculations  from  the  EPRI  Handbook.  The  agreement 
is  excellent 

For  the  experimental  work  described  below  a  constant  value  of  T)  was  used  for  all  tests. 
DE(T)  specimens  were  restricted  to  0.6  <  a/W  <  0.9  and  for  all  tests  T|  was  set  equal  to  0.27 
while  7  was  taken  as  (t|  - 1)  or  -0.73.  The  negative  y  did  not  have  a  strong  effect  on  these  tests 
because  of  the  small  amount  of  crack  extension  investigated  using  the  DE(T)  specimens. 

2.4.3  DE(T)  Crack  Length  Estimation 

The  DE(T)  specimen  is  tested  with  a  small  remaining  ligament,  generally  in  the  range  of 
0.6  <  a/W  <  0.9.  In  this  range  The  Handbook  compliance  equation  would  be  very  accurate  even 
for  the  pin  loaded  specimen  used  here.  The  compliance  equation  used  has  the  form; 

&/P  =  2^  (33) 

E' 


with: 


V(a/W)  =  —  {0.454sin«*-0.065sin’«*-0.007sin’u*+cosh->[sec(«')]}  (34) 

«* 


where: 


itfl 

2W 


(35) 


This  equation  must  be  inverted  to  be  used  for  unloading  compliance.  This  inversion  can  be 
performed  in  a  standard  fashion  to  give  a  polynomial  compliance  equation  of  the  form: 

a/W  =  0.0955026 -0.097503v'+0.245981(vf-0.115274(vy  (36) 

+0.0205763(vy  -0.0013593(vf 


with: 


(37) 
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Figure  6  Predicted  plastic  r\  factors  for  the  DE(T)  specimen. 
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2.5  Constraint  Correlations 


Recently  two  quantities  have  been  proposed  to  quantify  the  "constraint"  present  for  a 
given  combination  of  crack  geometry,  mode  of  loading,  and  material  toughness.  These  are  the 
'T-stress"  approach  (Bet6gon  and  Hancock  1992,  Al-Ani  and  Hancock  1991)  and  the  Q  parameter 
of  O’Dowd  and  Shih  (1991,1992).  These  two  quantities  are  described  and  utilized  separately  in 
the  next  two  subsections. 


2.5.1  T  Stress  Indexing  Parameter 

The  linear  elastic  crack  tip  stress  field  for  a  crack  along  the  negative  x-axis,  with  its  tip 
at  the  origin,  has  the  form: 


= 


fJB)  - 


(39) 


(40) 

V2itr 

where  the  T^  term  is  the  only  term  of  order  that  exists,  and  it  only  affects  a„.  Larsson  and 
Carlsson  (1973)  showed  that  the  sign  and  magnitude  of  this  term  does  alter  the  size  and  shape 
of  the  plastic  zone,  and  recently  Betegdn  and  Hancock  (1992)  and  Al-Ani  and  Hancock  (1991) 
have  suggested  that  the  amplitude  of  the  T,,  term  may  be  an  effective  constraint  indexing 
parameter,  even  in  the  elastic-plastic  regime.  They  show  that  low  constraint  geometries  like  short 
crack  or  tensilely  loaded  geometries  have  different  T„  values  than  deep  cracked  bend  geometries, 
and  they  suggest  that  the  T^  difference  causes  higher  apparent  toughnesses  to  be  found  in  such 
cases. 

To  assess  the  effect  of  T„,  a  biaxiality  parameter  is  used  having  the  form 

R  =  r  (41) 

This  quantity  has  been  evaluated  for  various  test  geometries  using  finite  element  and  other 
methods  and  is  available  in  the  literature  (Leevers  and  Radon  1982,  Kfouri  1986,  Sham  1991). 
Results  for  SE(B),  SE(T),  and  DE(T)  specimens  are  shown  in  Figure  7.  For  convenience,  this 
data  has  been  used  to  develop  polynomial  relationships  giving  P  in  terms  of  a/W: 
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SE(B): 


P  = 


SECT): 


P 


DE(D: 


p  =  -0.5844  +  0.6249 


-  1.4441 

(— f  +  11.264f- 
Wl  ( 

/ 

V  \ 

n\*  ( a\ 

-  20.950 

w]  ^ 

\ 

f  a  / 

-  1.6844 

■  1  * 

/ 

-  21.926[- 

V  "  /  \ 

-]  +  12.5641  f- 

)  -  1.352: 

1.2031 1 

7 

,  \3 


(42) 


(43) 


f-T 


(44) 


Using  these  relationships  the  value  of  ^  for  each  test  geometry  in  this  study  has  been 
calculated  as  shown  below.  Using  the  relationships  for  K  as  a  function  of  load  and  a/W  the  T,, 
value  at  the  crack  initiation  (J,J  point,  has  also  been  evaluated  for  each  specimen  tested.  These 
results  will  be  discussed  in  subsequent  sections. 

2.5.2  The  O  Indexing  Parameter 


The  Q  parameter  has  been  proposed  by  O’Dowd  and  Shih  (1991,1992)  as  an  extension 
of  the  Tg  concept  to  deformation  plasticity.  They  take  the  second  term  of  the  stress  fields  around 
the  crack  tip  in  a  power  law  hardening  material  and  propose  that  its  intensity  can  be  represented 
by  a  Q  parameter,  and  that  this  parameter  can  be  used  to  index  the  relative  constraint  of  a  test 
geometry.  The  authors  basically  relate  the  "true"  finite  body  stress  components  in  front  of  the 
crack  to  those  of  the  small  scale  yielding  (SSY),  infinite  body  case  in  the  form: 


=  ®  Jssr 

Q^o 

(45) 

*  ^jylssr  * 

(46) 

%  *  ®j^ssr 

(47) 

where  is  the  power  law  material  yield  stress  and  Q  is  a  dimensionless  parameter  found  by 
O’Dowd  and  Shih  to  be  between  0.2  and  -2,  in  front  of  the  crack.  Q  must  be  evaluated  using 
precise  finite  element  techniques  so  that  the  differences  between  the  true  stress  field  and  the  SSY 
stress  field  can  be  accurately  determined.  Results  for  Q  for  the  SE(B)  geometry  fi-om  O’Dowd 
and  Shih  (1992)  for  the  case  of  n  -  10  are  shown  in  Figure  8,  and  used  to  estimate  Q  values  for 


the  SE(B)  specimens  tested  in  this  program.  Q  values  for  the  SE(T)  and  DE(T)  specimens  have 
been  obtained  from  analysis  provided  by  Dodds'  as  shown  in  Figure  9-1 1  for  two  SE(T)  and  one 
DE(T)  geometry. 

2.6  Calculation  of  J.. 

In  order  to  compare  across  a  range  of  ^)ecimen  types,  special  care  had  to  be  taken  in 
its  evaluation.  The  1;^  value  calculated  is  dramatically  affected  by  the  initial  crack  length.  Ho, 
used  to  calculate  the  crack  extension,  Aa,  during  the  test  The  ASTM  Task  Group  E8.08.03  has 
proposed  a  method  to  utilize  the  initial  test  data  to  determine  a  best  a^  for  use  in  calculating  Aa 
and  brace  in  1,^  calculations^.  The  proposed  method  involves  fitting  a  straight  line  of  slope  ZOy 
to  the  initial  J  -  Aa  data,  and  choosing  the  best  fit  line  to  the  data  in  the  interval  0.2Jq  <  J  < 
0.6Jq.  Since  the  choice  of  a„  affects  the  subsequent  Jq  value,  an  iterative  process  is  necessary 
to  obtain  the  final  best  fit  Sq.  Jq,  and  hence  a  1,^  value  from  a  particular  experimental  data  set. 
The  proposed  method  also  adds  the  requirement  that  at  least  three  data  pairs  exist  in  the  region 
0.2Jq  <  j  <  0.6Jq  so  relatively  dense  data  is  required  for  a^  and  1,^  evaluation. 

This  method  was  initially  used  in  this  program,  but  it  was  found  that  the  method 
disqualified  some  specimens  due  to  a  lack  of  data  In  the  region  0.2Jq  <  J  <  0.6Jq.  In  place  of 
the  above  ASTM  method  an  alternative  method  was  developed  that  was  not  so  dependent  on  the 
early  data  on  the  J-R  curve.  In  this  method  the  relationship: 

flj  =  +  AJ/  + 

was  fit  to  the  {ij  -  Jj  data  of  each  set  from  the  minimum  a^  to  a^  +  2.5  mm  of  crack  growth,  as 
shown  in  Figure  12.  A  least  squares  procedure  was  used  to  evaluate  the  coefficients  a^.  A,  and 
B,  and  the  a^  parameter  was  the  desired,  best-fit,  initial  crack  length  which  was  then  used  to 
calculate  the  Aaj  quantities  for  the  J-R  curve.  This  method  worked  well  on  all  specimens 
analyzed  here  and  in  all  cases  adequate  data  was  available  for  the  fit,  no  iteration  was  needed, 
and  the  least  squares  technique  gave  a  unique  and  fully  defined  average  crack  length  for  each 
specimen,  while  the  proposed  ASTM  method  gives  only  a  range  from  which  each  investigator 
could  chose  a  somewhat  different  result  based  on  the  details  of  his  iteration  procedure. 


'Private  communication,  R.H.  Dodds,  Jr.  1993. 

^  "New  Standard  Method  for  J-Integral  Characterization  of  Fracture  Toughness,"  Draft 
Standard  of  ASTM  Subcommittee  E08.08,  March  1993,  American  Society  for  Testing  and 
Materials,  niiladelphia,  PA. 
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J/(a*Oo) 


Figure  9 


Q  =  (ayy-a,3y)/o„ 


The  Q  constraint  parameter  for  the  SE(T)  specimen  uiih  a/W=0.4,  nslO  (R.H. 
Dodds,  Private  Communication). 
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10  The  Q  constraint  parameter  for  the  SE(T)  specimen,  a/W=0.6,  n=10  (R.H.  Dodds, 
Private  communication). 
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Figure  11 
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-0.2 


-0.4 
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The  Q  constraint  parameter  for  the  DE(T)  specimen  with  a/Ws=0.7,  n=10  (R.H. 
Dodds,  Private  communication). 
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Figure  12  Polynomial  fit  of  crack  length  versus  J  for  determining  the  best-fit  initial  crack 
length,  a,. 
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3.1  The  Rotation  Correction 


The  SE(T)  specimens  were  a  new  geometry  and  problems  were  encountered  in  developing 
test  procedures  for  this  geometry.  One  obstinate  problem  encountered  in  testing  this  geometry 
was  the  initial  hook,  or  "crack  backup"  observed  in  the  J-R  curve  obtained  from  these  specimens 
as  shown  in  Figure  13.  The  unloading  compliance  measurements  show  a  clear  increase  in 
specimen  stiffness  as  one  brings  up  the  load  on  these  specimens,  both  in  the  elastic  and  the 
elastic-plastic  regimes.  This  problem  was  initially  assumed  to  result  from  the  round  hole  clevises 
that  were  used  for  the  first  tests  of  HY-100  (FYO).  New  clevises  were  machined  with  flat 
bottomed  holes  and  hardened  to  allow  a  free  rotation  of  the  specimen  halves  during  the  test.  This 
change  did  not  improve  the  results  however,  and  it  was  then  that  a  rotation  correction  was 
proposed  as  a  possible  solution.  The  rotation  correction  was  developed,  as  detailed  in  Section 
2.3.4,  in  a  fasldon  similar  to  that  used  in  ASTM  El  152  for  the  C(T)  specimen. 

The  effect  of  the  rotation  correction  was  to  remove  the  initial  crack  backup  without 
making  much  of  a  change  in  the  remainder  of  the  J-R  curve,  as  shown  in  two  examples  in 
Figure  14  and  Figure  15  The  most  dramatic  effea  is  at  the  start  of  the  J-R  curve  where  the 
apparent  initial  stiffening  of  the  specimen  is  corrected  to  show  an  improved  resistance  curve, 
especially  in  the  initial  portion,  which  improved  the  consistency  of  the  calculated  J,;  results. 

The  resistance  curve  slope,  T^,  measured  at  1  mm  of  crack  growth,  was  not  found  to  be 
changed  markedly  by  the  {^plication  of  the  rotation  correction,  being  beyond  where  most  of  the 
effects  of  the  rotation  correction  was  felt 

3.2  Ji-  and  T_..  Effects 

J,c  was  calculated  for  all  specimens  included  in  this  study  using  the  polynomial 
relationship  -  least  squares  fit  approach  described  above  to  determine  the  initial  crack  length,  Sq. 
The  results  for  both  materials  are  shown  in  Table  4  and  Table  5,  and  plotted  as  a  function  of 
specimen  type  in  Figure  16  and  Figure  17.  J,j  does  not  tqrpear,  from  these  plots,  to  be  very 
sensitive  to  the  specimen  type,  at  least  compand  to  the  large  variability  demonstrated  by  each 
specimen  type  taken  individually. 

The  material  tearing  resistance,  T^,  introduced  by  Paris,  et  al.  (1979),  was  also  evaluated 
for  these  specimens: 

T 

E  da 


where  Oy  is  the  material  yield  stress  and  dJ/da  is  evaluated  by  fitting  a  two  parameter  power  law, 
as  used  to  evaluated  Ji^,  to  the  J-R  curve  data  in  the  standa^  ASTM  E813  exclusion  zone,  and 
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kJ/m 


Crack  Extension  nun 


Figure  13  J-Resistance  curve  from  an  SE(T)  specimen  exhibiting  "crack  backup"  in  the  early 
portion  of  the  test. 
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Crack  Extension  mm 


Figure  14  Comparison  of  J-R  curves  for  an  HY-100  SE(T)  specimen  with  and  without 
rotation  correction. 


Figure  15  J-R  curve  for  an  HY-lOO  SE(T)  specimen  with  and  without  rotation  correction. 


Table  4 


Fracture  toughness  and  tearing  modulus  for  HY-100  specimens. 


Specimen 

ID 

Type 

a/W 

B 

(mm) 

4 

(kJ/mm^) 

T 

*  "ft 

(1  mm) 

FYO  1 

SE(B) 

0.66 

50. 

111.7 

25.2 

FY0  3 

SE(B) 

0.66 

50. 

137.8 

31.2 

FYO  21 

SE(B) 

0.14 

50. 

177.2 

49.9 

FYO  26 

SE(B) 

0.13 

25. 

163.7 

35.6 

FYO  27 

SE(B) 

0.14 

25. 

173.7 

48.0 

FYO  150 

SE(B) 

0.61 

25. 

145.1 

21.3 

FYO  151 

SE(B) 

0.61 

25. 

129.7 

26.6 

FYO  158 

SE(B) 

0.60 

12.5 

140.5 

23.8 

FYO  159 

SE(B) 

0.62 

12.5 

179.9 

24.2 

FYO  160 

SE(B) 

0.11 

12.5 

120.1 

38.2 

FYO  161 

SE(B) 

0.11 

12.5 

113.9 

40.9 

FY0  2SB 

SECT) 

0.40 

25. 

166.7 

47.9 

FY0  3SB 

SECT) 

0.47 

25. 

245.0 

50.1 

FY0  4SA 

SECT) 

0.65 

25. 

254.1 

31.5 

FYO  lOSA 

SECT) 

0.35 

25. 

201.8 

57.6 

FYO  USB 

DECT) 

0.68 

25. 

88.4 

31.0 

FYO  12SA 

DE(T) 

0.61 

25. 

107.1 

39.2 
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Table  5 


Fracture  toughness  and  tearing  modulus  for  A533B  specimens. 


Specimen  ID 

Type 

a/W 

4 

kJ/m* 

T 

^mat 

(1  mm) 

CT3 

C(T) 

0.6 

244.4 

122.4 

CT9 

C(T) 

0.6 

265.0 

125.3 

CTIO 

C(D 

0.6 

239.2 

112.1. 

DBl 

SE(B) 

0.62 

240.3 

149.9 

DB2 

SE(B) 

0.62 

308.5 

142.9 

DB3 

SE(B) 

0.62 

323.8 

160.5 

SBl 

SE(B) 

0.15 

306.3 

238.2 

SB2 

SE(B) 

0.15 

228.7 

260.8 

SB3 

SE(B) 

0.15 

333.0 

239.2 

SENl 

SECT) 

0.40 

145.5 

228.7 

SEN2 

SECT) 

0.36 

159.5 

289.8 

SEN4 

SECT) 

0.60 

305.2 

229.1 

SEN9 

SECT) 

0.40 

598.6 

288.4 

SENIO 

SECT) 

0.60 

231.8 

178.6 

SE5D 

DE(D 

0.7 

164.6 

301.1 

SE6D 

DE(T) 

0.7 

182.1 

257.3 

SE7D 

DE(T) 

0.68 

189.8 

253.3 

SE8D 

DE(T) 

0.71 

226.6 

255.0 

Jic  kJ/m 
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then  evaluating  dJ/da  at  a  crack  extension  of  1  nun  (0.04  in).  The  measured  values  obtained 
from  this  set  of  specimens  are  shown  in  Table  4  and  Table  5  and  plotted  as  a  function  of 
specimen  type  in  Figure  18  and  Figure  19.  The  tearing  resistance  varies  for  each  material  by  a 
factor  of  about  2.5  and  appears  to  be  dependent  on  the  specimen  type. 

In  the  following  sections  these  results  for  1,^  ^nd  T^  will  be  correlated  with  the  T^  and 
Q  constraint  parameters  introduced  in  Section  2. 

3.3  J-R  Curves 


Figure  20  and  Figure  21  show  the  measured  J-R  curves  for  the  baseline,  deep  notched 
SE(B)  and  C(T)  specimens  of  each  material.  The  variability  shown  in  these  figures  is  assumed 
to  be  due  to  material  variability,  and  is  typical  of  what  is  usually  found  for  structural  steels.  The 
dashed  bounding  lines  shown  in  the  figures  will  be  used  an  later  plots  for  comparison  of  the 
baseline  results  and  results  from  the  non-standard  specimens.  Comparisons  of  these  baseline  J-R 
curves  with  the  J-R  curves  of  the  short  cracked  SE(B),  SE(T),  and  DE(T)  specimens  for  both 
materials  are  shown  in  Figures  22-27.  The  most  immediate  observation  that  can  be  made  from 
these  figures  is  that  short  cracks  and  tensile  loading  seem  to  have  little  effect  on  the  Jj^  value, 
but  a  measurable  effect  on  the  slope  of  the  J-R  curve,  or  Tg,^  values,  with  higher  slopes  being 
found  for  all  of  the  short  crack  and  tensilely  loaded  specimens  in  con^arison  with  the  standard, 
deeply  notched  geometries. 

The  SE(T)  specimen,  SEN9,  appears  to  be  an  outlier,  with  crack  initiation  being  delayed 
for  some  unknown  reason.  This  specimen  has  been  explored  extensively,  but  no  reason  has  been 
found  for  its  elevated  toughness  behavior.  The  pre-crack  was  straight,  the  initial  and  final  crack 
lengths  were  estimated  accurately  by  the  compliance  technique,  and  the  load  displacement  curves 
seem  correct  in  every  respect  Figure  28  shows  a  plot  of  the  load  versus  COD  record  for 
specimen  SEN9  and  also  specimen  SENl,  which  had  a  nearly  identical  crack  length.  Specimen 
SEN9  is  clearly  much  tougher,  showing  a  continually  rising  load  displacement  record  throughout 
the  test,  while  specimen  SENl  rises  to  a  maximum  load  and  falls.  The  compliance  method 
predicted  crack  initiation  as  shown  on  Figure  28,  with  the  crack  initiation  of  specimen  SENl 
occurring  at  about  the  separation  point  of  the  two  curves,  while  the  crack  initiation  point  tor 
specimen  SEN9  was  much  later.  In  both  cases  the  extent  of  ductile  crack  growth  agreed  well 
with  the  post  test  optical  measurement 

3.4  Constraint  Correlations 

3.4.1  T,  Correlation 

The  Tg  quantity  was  calculated  for  each  specimen  from  the  tqiplicable  K,  a/W,  and  ^  at 
J,f  and  is  tabulated  in  Table  6  and  Table  7.  Figure  29  and  Figure  30  show  plots  of  J,c  versus  T^ 
for  each  material.  The  short  and  deep  SE(6)  specimens  for  both  materials,  and  the  C(T) 
specimens  for  the  A533B  alloy,  appear  to  be  insensitive  to  the  ^plied  T^.  The  tension-loaded 
SE(T)  and  DE(T)  specimens  are  in  general  agreement  with  the  exception  of  specimen  SEN9  for 
the  A533B  material,  which  was  discus^  previously  as  an  {q)parent  outlier.  The  HYIOO  tension- 
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Figure  20  Baseline  J-R  curves  for  HY-100  from  deep-cracked  SE(B)  specimens. 
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Figure  21  Baseline  J-R  curves  for  A533B  from  deep-cracked  SE(B)  and  CCD  specimens. 
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Figure  22  J-R  curves  for  HY-100  from  shon-cracked  SE(B)  specimens. 
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Figure  24  J-R  curves  for  HY-100  from  SE(T)  specimens  (with  rotation  correction). 
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Figure  25  J-R  curves  for  A533B  from  SE^  specimens  (with  rotation  conection). 


Figure  27  J.R  curves  for  A533B  from  DE(D  specimens. 
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A533B  SENS 


Load  versus  COD  curves  for  specimens  SEN9  and  SENl  of  A533B. 


Table  6 


Constraint  parameters  Tg,  and  Q  at  crack  initiation  in  the  HY-100  specimens. 


Specimen 

ID 


FYO  1 


FY0  3 


FYO  21 


Type 


SE(B) 


SE(B) 


SE(B) 


a/W 

B 

mm 

P 

T„@ 

MPa 

0.66 

50. 

0.476 

229. 

0.66 

50. 

0.476 

255. 

0.14 

50. 

-0.311 

-410. 

1  FYO  26 

SE(B) 

0.13 

25. 

-0.322 

-425. 

-0.8 

1  FYO  27 

SE(B) 

0.14 

25. 

-0.311 

-406. 

-0.8 

FYO  150 

SE(B) 

0.61 

25. 

0.395 

226. 

-0.1 

1  FYO  151 

SE(B) 

0.61 

25. 

0.395 

213. 

-0.1 

FYO  158 

SE(B) 

0.60 

125 

0.381 

216. 

-0.1 

FYO  159 

SE(B) 

0.62 

12.5 

0.411 

259. 

-0.1 

FYO  160 

SE(B) 

0.11 

12.5 

-0.345 

-423. 

-0.8 

FYO  161 

SE(B) 

0.11 

12.5 

-0.345 

-412. 

-0.8 

FY0  2SB 

SECT) 

0.40 

25. 

-0.270 

-182. 

-0.25 

FY0  3SB 

SECT) 

0.47 

25. 

-0.183 

-138 

-0.30 

FY0  4SA 

SECT) 

0.65 

25. 

-0.101 

-65.9 

-0.11 

FYO  lOSA 

SECT) 

0.35 

25. 

-0.325 

-259. 

-0.31 

FYO  USB 

DE(D 

0.68 

25. 

-0.407 

-217. 

-0.43 

FYO  12SA 

DEa) 

0.61 

25. 

-0.433 

-268. 

-0.40 

5 


Table  7 


Constraint  parameters  P,  T^,  and  Q  at  crack  initiation  for  the  A533B  specimens. 


1 

Specimen  ID 

Type 

a/W 

P 

To  @  J.C 
MPa 

CT3 

C(T) 

0.6 

0.573 

428. 

-0.1 

CT9 

C(T) 

0.6 

0.573 

446. 

-0.1 

CTIO 

C(T) 

0.6 

0.573 

423. 

-0.1 

DBl  SE(B)  0.62  0.411  300.  -0.1 


DB2 

SE(B) 

0.62 

0.411 

339. 

-0.1 

DB3 

SE(B) 

0.62 

0.411 

348. 

-0.1 

SBl 

SE(B) 

0.15 

-0.299 

-498. 

-0.8 

SB2 

SE(B) 

0.15 

-0.299 

-431. 

-0.8 

SB3 

SE(B) 

0.15 

-0.299 

-520. 

-0.8 

SENl 

SECT) 

0.40 

-0.271 

-172. 

-0.30 

SEN2 

SECT) 

0.36 

-0.315 

-220. 

-0.30 

SEN4 

SECT) 

0.60 

0.009 

6.7 

-0.22 

SEN9 

SE(T) 

0.40 

-0.271 

-348. 

-0.60 

SENIO 

SECT) 

0.60 

0.009 

5.8 

-0.20 

SE5D 

DE(D 

0.7 

-.0397 

-285. 

-0.70 

SE6D 

DE(T) 

0.7 

-0.397 

-300. 

-0.71 

SE7D 

DE(T) 

0.68 

-0.407 

-311. 

-0.72 

SE8D 

DE(T) 

0.71 

-0.392 

-335. 

-0.76 

Jic  kJ/m 
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Figure  29  Fracture  toughness.  as  a  function  of  T,  for  the  HY-100  specimens. 
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loaded  specimens  show  considerable  scatter,  with  the  SE(T)  specimens  being  high  relative  to  the 
SE(B)  results  and  the  DE(T)  specimens  being  low.  It  is  quite  possible  that  improvements  in  test 
technique  that  have  been  incorporated  in  the  more  recent  A533B  tests  in  part  explain  the 
reduction  in  scatter  shown  by  the  A533B  tensile  results,  in  comparison  with  the  earlier  HYIOO 
results. 


Hgure  31  and  Figure  32  show  plots  of  versus  Tg  for  each  material,  and  now  a  clear 
trend  of  material  tearing  resistance  versus  constraint  is  {q)parent  in  the  data,  with  the  tearing 
resistance  being  more  than  doubled  when  measured  with  the  low  constraint  short  crack  SE(B) 
specimen,  in  comparison  to  the  standard,  high  constraint,  deeply  notched  bend  or  compact 
specimens.  The  HYIOO  again  shows  much  more  scatter  in  values  than  does  the  A533B 
material.  This  might  be  due,  in  part,  to  the  much  smaller  values  obtained  from  the  HYIOO 
material,  as  well  as  from  basic  improvements  in  the  SE(T)  and  DE(T)  test  methods  between  the 
HYIOO  tests  and  the  AS33B  tests.  The  appearance  of  these  results  is  improved  by  the  lower 
scatter  for  each  specimen  type  demonstrated  by  the  Tg^  parameter,  in  comparison  with  the  1,^ 
parameter. 

3.4.2  O  Correlation 

The  Q  quantity  was  evaluated  for  each  specimen  from  the  applicable  analysis  using  Figs. 
8-10,  and  the  Ji^  and  a/W  and  is  tabulated  in  Table  6  and  Table  7  for  each  specimen  tested. 
Figure  33  and  Figure  34  show  plots  of  versus  Q  at  1,^  for  all  specimens  of  each  material.  As 
above,  neither  material  shows  dependence  of  on  Q,  at  least,  any  that  can  be  separated  from 
the  material  and  test  variability.  The  HYIOO  shows  much  more  scatter  than  is  shown  by  the 
A533B  material.  Overall  it  seems  that  is  not  strongly  affected  by  constraint,  as  measured  by 
Q,  in  these  specimen  geometries. 

Figure  35  and  Figure  36  show  plots  of  T,„„  versus  Q  (at  J^)  for  each  material.  A  clear 
trend  is  now  shown,  with  higher  constraint  resulting  in  lower  tearing  resistance  in  these  materials. 
For  both  materials  there  appears  to  be  a  rapid  increase  in  T,^  with  decreasing  Q,  then  a  leveling 
off  at  lower  Q’s. 
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Figure  31  Tearing  modulus,  as  a  function  of  for  the  HY-100  specimens. 
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Figure  32  Tearing  modulus,  T^,  as  a  function  of  T„  for  the  A533B  specimens. 
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Figure  33  Fracture  toughness,  J,.,  as  a  function  of  Q  for  the  HY-100  specimens. 
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Figure  34  Fracture  toughness,  as  a  function  of  Q  for  the  A533B  specimens. 
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Figure  35  Tearing  modulus,  as  a  function  of  Q  for  the  HY-100  specimens. 
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Figure  36  Tearing  modulus,  T„,  as  a  function  of  Q  for  the  A533B  specimens. 
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4.0  SUMMARY 


For  these  materials,  at  least,  the  toughness  parameter  is  relatively  insensitive  to 
constraint  as  characterized  by  or  Q,  while  the  tearing  resistance  is  elevated  by  reduced 
constraint  as  measured  by  either  or  Q.  Thus  the  J-Q  fracture  locus  is  inconsequential  for  fully 
ductile  behavior,  but  the  Tg,,,  -  Q  locus  is  an  important  curve  to  measure  the  extent  of  toughness 
enhancement  that  would  be  expected  under  conditions  of  low  constraint. 

The  best  results  for  a  low  constraint  specimen  geometry  are  obtained  from  the  short  crack 
bend  specimen.  This  specimen  is  now  relatively  easy  to  test  in  the  laboratory,  and  has  a  low 
constraint  by  any  measure,  especially  if  the  a/W  is  less  than  0.15.  Taken  together,  the  deep  and 
short  crack  bend  specimens  encompass  the  range  of  constraints  available  even  with  tension 
geomtries. 

It  has  been  demonstrated  (Kirk  and  Dodds,  1993)  that  the  plastic  ti  for  the  sort  crack  bend 
geometry  is  dependent  on  the  material  strain  hardening,  and  the  suggestion  has  been  made  in  the 
same  work  that  the  J  integral  should  be  measured  using  the  crack  mouth  opening  displacement 
using  Eq.  (9)  presented  earlier.  This  analysis  was  applied  to  the  specimens  of  this  program,  with 
a  typical  result  shown  in  Figure  37.  Three  methods  are  plotted  in  Figure  37.  The  first  analysis 
is  the  COD  J  integral  calculated  using  Eq.  (9),  the  second  analysis  uses  the  Sumpter  ti  analysis 
(1987),  but  does  not  apply  the  crack  growth  correction,  and  the  third  analysis  uses  the  Sumpter 
Ti  and  applies  the  crack  growth  correction,  the  method  used  in  the  work  presented  above.  In  the 
initiation  area  there  is  essentially  no  difference  between  the  three  methods.  Beyond  the  early 
region,  the  COD  J  integral  is  elevated  above  the  load-line  J  calculations.  In  the  later  part  of  the 
resistance  curve  the  two  non-crack  growth  corrected  methods  cross  over,  staying  above  the  crack 
growth  corrected  analysis.  The  maximum  difference  in  J  values  at  any  given  crack  extension  is 
on  the  order  of  10%.  It  is  not  clear  that  there  is  any  overriding  reason  to  prefer  any  one  of  the 
three  methods,  though  the  use  of  a  crack  growth  correction  has  strong  analytical  support,  and  has 
been  used  throughout  this  report. 
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Figure  37  Coo^arison  of  J-Resistance  curves  calculated  using  three  di^erent  J  formulations. 


S.0  CONCLUSIONS 


The  following  conclusions  are  drawn  from  this  work: 

1)  A  rotation  correction  is  essential  for  obtaining  an  accurate  J,c  or  J-R  curve  from  a  SE(T) 
specimen  of  the  type  used  in  this  work.  The  rotation  correction  developed  here  seems  to 
greatly  in^rrove  the  iq>pearance  J-R  curve. 

2)  J,g  does  not  seem  to  be  dependent  on  constraint,  as  applied  in  this  study,  at  least  as 
characterized  by  T^  or  Q. 

3)  The  material  tearing  resistance,  as  characterized  by  T^,  is  strongly  affected  by  constraint, 
with  T^  increasing  rapidly  with  decreasing  constraint.  This  is  true  whether  constraint 
is  measured  with  T^  or  Q. 

4)  The  best  low  constraint  test  specimen  geometry  is  the  short  crack  bend  specimen. 
Techniques  have  been  developed  which  make  it  relatively  easy  to  prepare  and  to  test 
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$|iccinwn  FYO-2Sn  JRAflleoulpul  luing Common  Method equaliom 


Specimen  l-YO-  3SH  jr/V  ni«  output  using  Common  Method  equations 


S  I  .£ 

OS'- 


^oeeoee>c>c>oc)^o^ooc>oo^oo 


^lc)^c>o‘oc>c>^oc>c»c>^c»c>o^o^ooooddo 


9  B 

Si 

=  ®'«RRp:8Sa|KffS 

^  s 


S  Q  s  3;  s 
§  »  S  s  »  K 


-sl'rJisississoosisssssis&sssqs::; 

S  S.^6c>ee>e»c>oe>e>e»e>e>de>beeoeeoooooo 
U  5  III 


C  iS  - 


iTi'Sr-'r-lQ^ovp  —  — 


,  SKr^82t;saa 


I  .s 


S50  —  •o«r.  — 


1 

£j^^«SPSSgg||g| 

8  I  i( 

w 

iJs;ig||i§ssss 


«  SO  r«  r^  80  eo 


•■•^^^♦^«©r<8eaopo*c.  f*‘r>4 
—  flO«r.  ^®805p^^ 
ooOkOO«»«^«nv%>e'2r*^9o^o 


ste:?3  =  S!(-;(;ie£3:£S; 

r^d0^^c^^C^^^80^>»y5'«'H^  fn 


od^oddedo 


<ri  «r.  «r.  «r.  «  50  ^  r* 

•J  v?97!3i9S|ssS^|a 
|JS|s:::2SSSRSS82Z 


f^«%r<sfseo«on 

9:SR7S&8r:P.  9;zs;9?8! 

Q&®®080«or->C¥*.  »r.  f^fSfNr* 
«  «r)  IT,  «r.  •r.  irv  tr.  »r,  ir,  irv  m“,  ¥%  *r.  ir. 


-  l*^gggO  —  so 

'•S.  &  ‘^SS88ooeeooSS 
^  .:,  d  d  d  d  d  d  d  d  d  d  d  d 

5  c  ^ 


r*80«0^O--— 

oSeoSSSSSSSSSS 

d  d  d  d  d  d  d  d  d  d  d  d  d  d 


oc  r-  V  flc  r* 


Bl  5 

w  w 


lllllllllllll 

ddddddddedddd 


JS?!SrSgC:P8SS?; 

O  2r*5^»^«or»«r.  r*i«^d 

K  •»  £  M-.  »r.  ir,  ir.  IT.  IT.  ^  ^  ^  f»^  ^ 


ss;^7;^iesss289::‘ 

=  C9;Q||?g&8S§»o 

n  «r.  ^  qo  ^  fM  »r.  ^  ^  fs 

n•-ddpo»pcpr•r•«*^•r.  *w 


Os«r.  osir.  — 

«r.  ®fs|so<^  —  tr,r^ot  —  f^ 

c-^'ggSpoSSgSSSo 

0.S!dddddddddddd 


ir.  «r.  5*  «  IT.  o>  r>*  ^  C.  r«i  Ov  rq  IT. 

r«*OvO«Nf^«r.  so 
fn  ^  ^  ^  ^  ^  <r,  »r.  »r,  sr, 

ooooooeooooooo 

dddddddddddddd 


•8  2  li 

i^i\ 


SRs8isis:!:!^S7SS 


«r>mp«r.  r>i«09>r4^se^  —  cetr.  Q 
»^^8or>—  o««C80^<Fr«>«aod 


—  iNPT^ir.  ^r^®0^0  —  iNfn^ir.  sor*80O'C  — 


JKA  file  oulptii  using  Common  Method  equations 
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S|Kcinien  TYO-  KISA  JKA  Tile  nuipul  using  ComnHMi  Method  equations 
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Unknd  Na  Loadline  CX)D  Plastic  Lmd  Loadline  COD  Della  J  J 

No.  Data  Slope  Slope  Corr.  Area  a  Plastic 
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Specimen  CT-3  JRA  Tile  milpul  using  Common  Method  equations 
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Specimen  Cf-  9  jHA  Tile  oulpul  using  Common  Method  equations 
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Spe*'"**"  CT- 10  JRA  nie  oulpul  tiling  Common  Method  equitioni 
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Specimen  Cl  - 10  JRA  file  oulpul  luing  Common  Method  equalioM 
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Specimen  DB-2  JRA  file  output  using  Common  Method  equations 
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JRA  file  oulpul  wing  ComnHM  Method  equalione 
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